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We consider the response of etalons created by a combination of a conventional mirror and a metasurface 
composed of a periodic lattice of metal scatterers with a resonant response. This geometry has been used 
previously for hybridization of localized plasmons with Fabry-Perot resonances and for perfect absorption, in 
so-called Salisbury screens. The particular aspect we address is if one can assume an environment-independent 
reflectivity for the metasurface when calculating the reflectivity of the composite system, as in a standard Fabry- 
Perot analysis, or whether the fact that the metasurface interacts with its own mirror image renormalizes its 
response. Using lattice sum theory, we take into account all possible retarded dipole-dipole interactions of 
scatterers in the metasurface amongst each other, and through the mirror. We show that while a layer-by-layer 
Fabry-Perot formalism captures the main qualitative features of metasurface etalons, in fact the mirror modifles 
both the polarizability and reflectivity of the metasurface in a fashion that is akin to Drexhage’s modiflcation of 
the radiative properties of a single dipole. 


I. INTRODUCTION 

Local enhancement of electromagnetic fields by periodic 
lattices of metal scatterers has many a ppli cations, which in¬ 
clude improving refractive index sensors^, brightening solid- 
state light emitterJ®, and providing substrates for linear and 
nonlinear optical spectroscopieP. These applications on one 
hand employ the fact that single metal objects can provide 
strong scattering resonances on basis of plasmonic (i.e., free 
electron) oscillations. In addition, these o bjec ts can be arrayed 
either in a diffractive “plasmonic crystal’® or can be packed 
more densely in a “metasurface” or “frequency-selective” sur- 
facP. Several works have investigated how the response of 
such periodic plasmonic antenna arrays can be further mod¬ 
ified by placing them close to a “ground plane”, i.e., a re- 
fiective interface. In fact, this idea dates to the 1950’s when 
the “Salisbury screen” was proposed as a highly efficient mi¬ 
crowave absorber in radar applicationP. This absorber is cre¬ 
ated by placing a weakly reflective and slightly absorptive 
layer at a quarter wavelength from a mirror. The concept has 
recently been extended, for instance, to mid-IR absorption in 
graphenP, and, in an independent approach, to metamaterial 
arrays at < A/10 distances from ground planes to reach above 
97% absorption in the mid-IR and THz rang^^^^. Ameling 
and Giessen have also pursued plasmonic arrays near refiec- 
tive interfaces in an effort to enhance plasmonic sensors by 
strongly coupling localized plasmons with Fabry-Perot reso- 
nanceP®^. These approaches combine two essential ingredi¬ 
ents: (1) a resonant scattering layer that provides a reflection 
constant going through a 180° phase shift and an amplitude 
maximum, and (2) tuning of a geometrical resonance associ¬ 
ated with the distance to one, or two, refiectors. 

While the structures outlined above have been the subject of 
ample numerical simulations, we have found it instructive to 
re-analyze their physics from two different perspectives. First, 
we note that the geometry of an etalon with a weakly refiec- 
tive interface in front of a unit refiectivity interface is of inter¬ 
est not just as an absorber, but also for its phase properties. In 
particular, this geometry, known as a Gires-Toumois etalon, 
provides large chromatic dispersion even with a non-resonant 
front interfacJ^. Here we examine the phase response tak¬ 


ing a resonant metasurface as front interface. Second, in addi¬ 
tion to reporting on the phase behavior of metasurface etalons, 
we are interested in determining how strongly the metasurface 
scatterers become polarized, and what the role of backaction 
is. It has been previously assumed that it is possible to assign a 
dispersive refiection constant to the bare metasurface, to then 
apply a Fabry-Perot multilayer refiectivity formalism to pre¬ 
dict the refiectivity of a composite system. This assumption 
is difficult to reconcile with the established notion of back- 
action, whereby the scattering strength and linewidth of plas¬ 
monic scatterers can be modified by refiective interfaces due 
to interference of an induced dipole with its own mirror im- 
agjmn. 

We report on a self consistent lattice sum theory that in¬ 
cludes backaction, as well as all dipole-dip ole in teractions of 
scatterers in the metasurface, both directl}P3^ and through 
multiple scattering off the nearby refiective interface. This 
serves as benchmark for the Fabry-Perot multilayer refiectiv¬ 
ity model postulated by Ameling et alP. We find that the 
Fabry-Perot model indeed gives an excellent estimate of the 
composite etalon performance. In particular, for dense lat¬ 
tices (pitch < A/5) and large distances from lattice to inter¬ 
face (> A/5), the model gives quantitatively accurate results. 
For more dilute lattices, the Fabry-Perot model reproduces all 
qualitative features, but refiection phase and amplitude can be 
off by as much as ~ 0.5 rad and 0.1 respectively. These devia¬ 
tions show the effect of backaction, or equivalently, of multi¬ 
ple scattering interactions with the mirror via evanescent grat¬ 
ing orders. In Fabry-Perot terms, the equivalent metasurface 
refiection constant is not independent of its environment. This 
result indicates that the response of stratified metasurface sys¬ 
tems cannot be rigorously captured in an effective multilayer 
formalism, even if in particular practical cases the multilayer 
formalism is appropriate. 


II. FABRY-PEROT ANALYSIS OF METASURFACE 
ETALONS 

To assess the general performance of etalons with a mir¬ 
ror and a metasurface interface, we first set up a simple two- 
interface model, similar to the efforts by Ameling et 
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We consider incidence from medium 1, onto a metasurface 
that is separated by a distance d from a highly reflective in¬ 
terface (medium between the two reflective planes labeled 2, 
beyond that, medium 3, refractive indices denoted as n^.)- We 
denote with k = uj/c the free space wavenumber, and indicate 
with Trim 5 tnm the reflection and transmission amplitude co¬ 
efficients when going from medium n to m. The reflectivity 
of the stack, under the assumption that a simple Fabry-Perot 
model is applicable, given by 


FP _ t2ir23h2 

•stack -ri2+ 


( 1 ) 


where r 23 is the reflectivity of the back reflector, taken from 
Fresnel equations for the reflectivity of a real metal mirror, 
while ri 2 , j 2 i and ti 2 , ^21 refer to the metasurface. Garcia de 
AbajcP^has worked out the reflectivity of a periodic surface 
of scatterers in free space, under normal incidence as r = ri 2 = 
^ 21 , t = ti2 = ^ 21 , 


27rik 1 

(2^ 1 /Clfree “ Q 


with t = 1 + r (2) 


with afree the particle polarizability (for a single scatterer, in 
absence of the interface), the area of the unit cell, and Q 
an Ewald lattice sum that takes into account all dipole inter¬ 
actions. We refer to Ref. [24] and the Appendix]^ for the def¬ 
inition and evaluation of free space lattice sums. Let us for 
simplicity assume plasmon nanobars, that is, particles that are 
only polarizable along one of their axis (x as the incident po¬ 
larization) and along one axis of the lattice, which applies for 
both the experiments of Ameling et al.^, and the experiment 
described in Appendix[^(performed with split rings with neg¬ 
ligible magnetic moment and for which, at normal incidence, 
the magnetic dipole does not contribute to reflection and scat¬ 
tering because is out of plane). According to Eqs. (10) and 
(13) of Ref. |24l for cubic, nondiffractive lattices with X > a 
the sum is well approximated by 


^ / 47r^V2 1 

\ y27r/(fca) 

If the electrostatic polarizability ao is chosen LorentziarP^I^, 
the full particle polarizability including radiation damping 
(also known as ‘dynamic polarizability’) reads 



ao = 


Vco^o 


cJq - i^uj -i^k^ 


( 4 ) 


with the term containing k = nujic representing radiation 
loss. Thus an explicit expression results for a resonant re¬ 
flectivity around ccq with Ohmic damping rate 7, and with a 
strength tuned through scatterer volume V. We discuss re¬ 
sults for a typical scenario with scatterers resonant at tele¬ 
com wavelengths (1.25 x 10^^ s“^), taking 7 = O.Iccq typ¬ 
ical of the Ohmic damping in gold. We choose a volume 
V = 0.0020 jjivc? to obtain an extinction cross section on res¬ 
onance {AiTklma) halfway the maximum cross section for a 
strong dipole scatterei^(3/27rA^), as is typical for plasmonic 
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FIG. 1. Amplitude and phase of the stack reflectivity calculated using 
a simple Fabry-Perot model for dilute to dense systems. The pitches 
are indicated in the figure and the respective unit cell areas are (a) 
o? = 1 fim^, (b) o? = 0.25 jim^ and (c) = 0.09 jim?. 


scatterer according to measurements in Ref. |28| To quantify 
r 23 we take the Fresnel coefficients for an air-silver interface, 
using tabulated optical constants for silvei^S. We note that 
this Fabry-Perot model only uses normal-incidence reflectivi¬ 
ties of the interfaces, and neglects possible grating diffraction 
effects. While in this paper we consider structures that have no 
far-field diffraction orders in the frequency range of interest, 
the full model discussed in section illll includes all diffraction 
orders. 

Figure shows the expected amplitude and phase response 
for three cases ranging from a dense to a quite dilute lat¬ 
tice, where the bare lattice reflectivity goes from being much 
broader than 7 and strong (|rp > 80%, with < 1% transmis¬ 
sion and 20% absorption), to being limited in width to 7 and 
weak (|rp ~ 3% at ~ 70%). The etalon behavior is very 
different for these three cases. In the very dense case, the 
reflectivity shows very sharp bands of absorption around the 
integer roundtrip condition, with a lineshape that is asymmet¬ 
ric, and reverses asymmetry when uj sweeps through ujq. The 
very narrow slivers of near-perfect absorption in uj - d space 
are associated with very sudden changes in phase of the re¬ 
flectivity. As the lattice is diluted, the narrow stripes of ab¬ 
sorption broaden in d-range and narrow in cc-range. Points 
of perfect absorption occur on either side of ujq, at opposite 
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sides of the geometrical resonance condition, and are associ¬ 
ated with phase singularities in the reflectivity. Finally, for the 
very dilute case, as function of spacing to the mirror pock¬ 
ets of near-perfect reflectivity and strong absorption alternate, 
while no phase singularities occur. The following remarkable 
features stand out. First, perfect reflectivity can be obtained 
even though the single pass absorption by the metasurface is 
between 20% and 50% for the dense to dilute lattices. Sec¬ 
ond, for the dense lattice the very broad and quite absorptive 
bare response for the lattice results in very sharp asymmetric 
response line^. Third, for the dilute lattices very strong re¬ 
flection contrast is possible, even though the metasurface itself 
only has a few percent reflection. In part this is understand¬ 
able from the fact that the single pass absorption is quite high 
(30% for the most dilute lattice). This k survey shows that 
a Fabry-Perot analysis of metasurface etalon behavior con¬ 
tains the physics of the Salisbury screen (highly reduced re¬ 
flection and increased absorption), the Gires-Tournois etalon 
(strong chromatic dispersion even in wavelength ranges where 
ri 2 does not change phase), and the apparent transition from 
crossing to anticrossing behavior of the Fabry-Perot disper sive 
resonance and the LSPR reported by Ameling et al.^^ElISl 


III. FULL THEORY 


Having introduced the physics of metasurface etalons in a 
Fabry-Perot model, we now turn to a full point-dipole model 
to confirm when the accepted Fabry-Perot intuition is correct, 
and to identify the origin of any deviations from it. Our goal is 
to And the full lattice-summed Green function Q for a system 
composed of a infinite two-dimensional lattice of point scat- 
terers at a distance d from a parallel interface, taking into ac¬ 
count all the multiple scattering interactions among particles, 
including those via the mirror. While we only write down the 
theory for electric dipole moments p, the model directly ex¬ 
tends to magnetic and magnetoelectric dipoles (read (p, m) 
for p in all equations, driven by (E, H) instead of just E, us¬ 
ing the 6 X 6 polarizabilites and Green functions of Ref.|23]and 
[30b . As in Ref. [30l our starting point is that the response of 
a periodic lattice of polarizable point particles illuminated by 
a plane wave (electric held Ee^^H '^ with parallel wave vector 
k||, plus its Fresnel reflection at the mirror) can be summa¬ 
rized entirely by the response p of the parti cle at the lattice 
origin (dynamic polarizability a) which readJ^^^^ 


p = cx 


E + ^G((Rn,d),(0,d))e'^i|-^-p 

ntO 


(5) 


so that 




G((R„, d), (0, d))e 


ikii-R, 


-E. (6) 


Here, the real space lattice vectors are taken parallel to the 
(x,^)-plane, as well as the interface which is at z = 0. With 
G(r,r') we denote the dyadic Green function in absence of 
the lattice, but accounting for the interface, for a single dipole, 
explicitly defined in Appendix]^ The polarizability a must be 


corrected for the presence of the substrate. This means it must 
be chosen to represent a self-consistent scatterer in presence 
of the environment specified by G(r, r'). Mathematically this 
means that the polarizability must be corrected for backaction 
through 


a ^ =afrgg-G((0,d),(0,d)) 

where is the dynamic polarizability in free space and 
G((0, d), (0, d)) represents the interface-induced backaction. 
The dyadic Green function for the case of a single planar in¬ 
terface can be separated into a free space and reflected parP3 

G(r,r') = Gfree(A^') + Grefl(r,r'). (7) 

The free space part is simply the held of a dipole in free space 
(same index as the medium in which the lattice sits), while 
Grefl is a more complicated function constructed in the an¬ 
gular spectrum representation, and requiring a parallel wave 
vector integral including the Fresnel coefficients of the inter¬ 
face system, described in detail in Appendix This integral 
has poor convergence properties that req uire complex integra¬ 
tion methods outlined by Paulus et al.^ 

Required for evaluating the lattice response through Eq. 
is the lattice-summed Green function defined through 

ntO 

where the superscript ^ denotes a sum over all terms n ^ 0, 
while lattice sums that do include the n ^ 0 term are denoted 
without superscript in this paper. Since the Green function in 
presence of the single planar interface separates as a sum, also 
the lattice-summed Green function can be written as a sum of 
a free space and a reflected part 

a*(k||,r) = C?;,33(k||,r) + arefl(k||,r). (8) 

The lattice summation req uired to obtain Qf^ee was discussed 
in detail in previous reportRegarding the reflected part, 
we And that there is an intuitive interpretation in terms of mir¬ 
ror dipoles. If one considers that the interface-corrected po¬ 
larizability is set by ol~^ = - G((0, d), (0, d)), one can 

rewrite Eq. as 

P = —^^ -E (9) 

<e-^(k||.r) 


with 


g(k||,r) = afree(k||,r) + erefl(k||,r). (10) 

In terms of mirror image dipoles, the interpretation is that a 
given particle interacts with all the particles except itself in the 
physical lattice (captured in and furthermore with all 

the mirror imaged scatterers including its own mirror image 
(captured in ^refi)- 

As regards rigorous evaluation of the reflected lattice sum, 

0refl(k||,r) = y) Grefl (Rmn-r II (H) 

m,n 
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the problem looks quite formidable since the reflected Green 
function Grefl(i*, r') that one needs to sum has the following 
mathematical structure, at distance d from the interface 


1000 nm pitch 500 nm pitch 300 nm pitch 



Here Mrefi depends on in-plane wavevector q and contains 
the Fresnel coefficients of the interfac^^, specified in Ap¬ 
pendix The Sommerfeld integral is subtle to perform, as 
a pole occurs at the light line due to the 1/kz term (where 
kz = - |qp represents out-of-plane wave vector), as well 

as at any wave vector q that corresponds to a guided mode at 
the interface. While summing an infinite set of Sommerfeld 
integrals appears daunting, in fact the completeness relation 

( 12 ) 

m,n ^ m,n 



leads to a remarkable simplification. Indeed substituting the 
completeness relation removes the Sommerfeld integration al¬ 
together, instead resulting in a quickly converging sum. Thus, 
owing to parallel wave vector conservation, the interface lat¬ 
tice sum is easier to calculate than the interface Green func¬ 
tion Grefl itself. The sum reads (with kf - ^Jk^ - |k|| + gp) 



distance (ijm) distance (ijm) distance (ijm) 


erefl(k||,r) = ^ E (^3^ 

^ m,n 

This sum runs over propagating as well as evanescent diffrac¬ 
tion orders, meaning that also any diffractive coupling to 
guided modes of the underlying mirror are accounted for, if 
relevant. 

Far-field observables, i.e., reflection, transmission and 
diffraction efficiencies follow simply once the induced polar¬ 
izability has been solved for, by using the far field asymptotic 
expansion of the reflected part of the Green function (Chapter 
10 in Ref. [ST]). As in Ref.[30l the far-field response consists 
of the specularly reflected beam as well as diffracted orders in 
the upper and the lower medium. We note that none of this 
analysis was specific to a substrate performing as a good mir¬ 
ror. In fact, this formalism can deal with lattices on arbitrary 
metallic or dielectric stratified media, where the lattices can 
be either diffractive or subdiffractive. The model also directly 
applies to, for instance, analyzing the strength of diffractive 
surface lattice resonances of lattices near the interface of two 
transparent media, as analyzed by Auguie et al.l^. Arbitrary 
stratified systems can be included by replacing the Fresnel co¬ 
efficients implicitly appearing in the reflected Green function 
by appropriate multilayer Fresnel coefficients. 

IV. RESULTS OF THE FULL THEORY 
A. Reflectivity 

Here we discuss reflection coefficients calculated through 
the exact Green function formalism set up in the previous sec¬ 
tion. For comparison. Fig. shows the absolute squared of 


FIG. 2. The reflection coefflcient for lattices of plasmon bars 
(1000,500,300 nm pitch going from left to right) from full lattice- 
sum point-dipole calculations, (a) Lattice reflectivity in absence 
of any mirror (thereby independent of distance), (b) full reflectiv¬ 
ity of the system calculated in presence of the mirror, (c) approxi¬ 
mate etalon reflectivity calculated by inserting the result of (a) in the 
Fabry-Perot equation, and (d) apparent lattice reflectivity retrieved by 
inverting the Fabry-Perot equation taking as input the mirror Fresnel 
coefflcient and the full reflectivity shown in (b). 


the frequency dependent reflection coefficient from point 
dipole (PD) lattice-sum calculations in absence of the mirror, 
Fig.|^(a), and calculations that contain the mirror, Fig.[^(b), 
as in the preliminary Fabry-Perot analysis. 

For relatively dilute lattices, e.g., pitch 500 nm, the abso¬ 
lute squared value of the reflection amplitude is small 

in the whole frequency range. Fig. (a). For the full sys¬ 
tem including mirror, the off-resonant |^i™ticeP almost 
100% throughout the distance range as expected since at 
these frequencies the response is dominated by the reflec¬ 
tion of the silver mirror. Fig. (b). On-resonance, we find 
a low reflectivity of less than 10% in minima at distinct 
mirror-metasurface distances. The combined system has a re¬ 
markably large reflectivity contrast for the reported frequency 
range of [1%,95%], even though the individual interfaces 
both feature a small reflectivity contrast of [0.5%, 5.5%] for 
the lattice and [98%, 99%] for the mirror. The fact that fringes 
occur at A/2 intervals is consistent with the idea of a Fabry- 
Perot resonator governed by constructive interferences when 
the round trip accumulates a full wavelength (2kd = m27r). 
At all densities, a fringe of perfect reflectivity occurs when 
the lattice is right in the node of the standing wave that forms 
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at the mirror. 

For more dense lattices, the etalon reflectivity transits into 
sharp dispersive features in the reflection coefficient. For the 
lattice with a pitch of 300 nm, the reflectivity is qualitatively 
very different, having gone from absorptive features on a high 
background to asymmetric dispersive lineshapes. Fig. |2](b) 
right column. The reversal of contrast at 1.2 x 10^^rads 
i.e., whether the reflectivity is low or high to the left, or the 
right of the Fabry-Perot resonance is due to the lattice reso¬ 
nance, whose response function goes through a tt phase flip. 
This transition could be understood as a transition from coher¬ 
ent mixing of a continuum (mirror reflectivity) with a weak 
absorption resonance (dilute lattices, walking along paths of 
constant separation), to mixing of a continuum (mirror reflec¬ 
tivity) with a strong resonance (lattice response). Thus the 
lineshape change is akin to that in the families of lineshapes 
of the Fano interference effect that has attracted signiflcant 
attention in marw communities, including recently that of op¬ 
tical scattering^. 


B. Fabry-Perot model and backaction effect on reflectivity 


We And that the Fabry-Perot model is remarkably effective 
for predicting the reflectivity of the metasurface etalon on ba¬ 
sis of just the mirror reflectivity, and the metasurface normal- 
incidence reflectivity. To bring this out, we compare the full 
model (i.e., including all particle-particle interactions, includ¬ 
ing those through the interface), with the Fabry-Perot model 
that takes as only input the mirror reflectivity (in absence of 
the lattice), and the bare lattice reflectivity (i.e., reflectivity 
calculated with the point dipole method in absence of the sub¬ 
strate). The resulting is shown as panel (c) for each 

lattice density in Fig. The same three regimes of reflec¬ 
tivity depending on the strength of the metasurface response 
are seen, with very similar magnitudes of reflectivity features. 
Inspection shows that the calculated reflectivity features are 
qualitatively very similar. An important result of our work is 
that it verifles that the Fabry-Perot model is highly appropriate 
for metasurface etalons. 

More thorough analysis in fact does show deviations. For 
instance, the precise location of troughs in the reflectivity 
shifts as much as 100 nm in distance and up to 5% in fre¬ 
quency, with differences of order 0.1 in reflectivity contrast. 
A clearer representation is offered by constant-frequency line- 
cuts that show the reflection amplitude, and the reflection 
phase in both models. Figureshows such cuts at cc = ccq. At 
pitches of 500 (A/3) and 1000 nm (2A/3) (not diffractive but 
not deeply subwavelength) the reflection amplitude differs by 
up to 0.1, while the reflection phase differs by up to 0.5 radi¬ 
ans (about 30°). This would constitute a signiflcant difference 
when using the stack as a phase-control element. Comparing 
line cuts at various frequencies, it is clear that the agreement 
is worst (i) with increasing pitch, being noticeable for pitches 
> A/5, (ii) closer to the mirror, and (iii) in the low-reflectivity 
areas, as expected given Eq. (IT ). Deviations from the Fabry- 
Perot model are due to evanescent diffractive contributions, 
which require signiflcant favoring small mirror-lattice 


d=1000 nm d=500 d=300 



FIG. 3. Reflection amplitude |r|and phase A(/? versus mirror- 
metasurface spacing for different lattice densities, at frequency uj = 
ujQ. Blue curves show the full point dipole result, while red curves 
show the Fabry-Perot result. It should be noted that to focus on the 
differences between the model, the plotted phase is relative to a sub¬ 
tracted trivial linear behavior (due to the fact that in scanning dis¬ 
tance, the front reflective surface approaches the observer). 


separation d, and large pitch (reducing ikf). At low reflectiv¬ 
ity the lattice is most strongly polarized, as at the lines reach 
unit efficiency, the lattice is at a held antinode, at which point 
(evanescent) diffractive corrections are irrelevant. The reason¬ 
ing in terms of evanescent diffractive orders is similar to that 
by Zhao et al.^S, where it was found that a similar criterion 
determines when stacked metamaterial layers are sensitive, or 
immune to lateral alignment errors. 

One viewpoint on the deviations between Fabry-Perot 
model and the model including all particle and mirror inter¬ 
actions, is that the metasurface reflectivity is apparently not 
independent of environment. A fundamental tenet of the stan¬ 
dard treatment of dielectric multilayers is that the composite 
reflectivity can be expressed rigorously in terms of just Fres¬ 
nel coefficients for reflection and transmission at each inter¬ 
face and phase terms associated with propagation across the 
thickness of each layer. As example, in the two-layer formula 
Eq. 0, the coefficients (r^m^^nm) are independent of the 
spacing d. The coefficients are determined for each boundary 
by boundary condition matching. Since a metasurface is not 
a microscopically smooth interface to which boundary con¬ 
dition matching applies the assumption that ri ,2 is indepen¬ 
dent of the substrate might not hold. We try then to determine 
what is the apparent lattice reflectivity that one needs to insert 
in the Fabry-Perot model to make it match the full calcula¬ 
tion. This is motivated by the use of retrieval procedures in 
metamaterial research, that aim to retrieve effective medium 
parameters by inverting experimental observations. Similarly, 
given the mirror reflectivity, and the mirror-to-lattice spacing 
{d and r 2 , 3 ), one can envision a retrieval procedure to deter¬ 
mine the apparent metasurface reflectivity from the measured 
complex etalon-reflectivity, which can be measures using an 
interferometer. Equation can be inverted by simple alge¬ 
bra to obtain a complex effective ri ,2 with no ambiguity for 
each separation d. Figure panel (d) shows the apparent lat¬ 
tice reflectivity for each lattice density. Clearly the retrieval 
is successful, in that the retrieved apparent reflectivity does 
not show the very strong fringing contrast of the composite. 
Nonetheless the apparent reflectivity is not identical to the re¬ 
flectivity of the lattice in absence of the mirror. Instead it does 



















6 


depend on the spacing d between mirror and lattice. The reso¬ 
nance frequency oscillates with distance to the mirror, as does 
the width. 

There are two ways to understand the dependance of the ef¬ 
fective reflectivity of the lattice on distance to the interface. In 
grating terminology, slicing a structure up in layers is only a 
valid approach if one includes sufficient evanescent diffracted 
orders to capture near field coupling. It is evident that a Fabry- 
Perot model that takes only the refiectivity and transmission 
of propagating orders does not include the evanescent orders, 
making it fundamentally impossible to predict the response 
of metasurface stacks using a simple Fabry-Perot model. As 
a second way to understand the variation of the apparent lat¬ 
tice refiectivity with the external geometry, we propose that 
it is akin to backaction correction to the response of a scat¬ 
tered Since the experiments by Drexhage in the 1960’s, it 
is well known that the radiative decay rate of quantum emit¬ 
ters, as well as the radiative damping rate of classical dipoles, 
depends on distance to a mirror, because interactions of the 
dipole with its own mirror image cause super- or subradiant 
dampin g. I t has been realized since an experiment by Buch- 
ler et al.ES that for single plasmon particles this means that 
the scattering cross section varies in strength, width, and reso¬ 
nance center frequency as the particle is held in front of a mir¬ 
ror (oscillation period around A/2). Here a similar back action 
effect occurs, however now in form of a lattice of dipoles in¬ 
teracting with an entire lattice of image dipoles. 


C. Backaction in polarizability 
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Besides analyzing back action effects in the observable, i.e., 
refiectivity, it is possible to determine how strongly the polar¬ 
izability of particles is renormalized. Underlying the varia¬ 
tion in scattering cross section of a single particle as it is held 
in front of a refiective surface, is the notion that its polariz¬ 
ability a is renormalized by the interaction with the mirror. 
The polarizability quantifies the strength of the induced dipole 
moment given that the incident field applies a set strength of 
1 V/m. In other words, if one has a single particle in front of 
a mirror, then upon driving by an external field Eq (in itself 
satisfying Maxwell’s equations in presence of the mirror, but 
absence of the particle), the expected dipole response is 

p = q;[Eo + Grefl(r,r)p] 

If one defines the polarizability as the response to the exter¬ 
nally applied driving field Eq, rather then the total field, the 
consequence is that the polarizability is renormalized 

p = O^effEo 

with 

— = — -Grefl(r,r). 

^eff ^free 

In a similar way, Eq. 0 shows that the per-particle polar¬ 
izability — when defined as the induced dipole strength in 


FIG. 4. Point dipole lattice-sum calculations to retrieve the polar¬ 
izability of nanorod scatterers at a distance from a silver interface 
driven by an x-polarized plane wave under normal incidence. The 
single-particle polarizability (a) varies with distance due to the in¬ 
teraction with its mirror image. Lattice arrangements of 1000 nm 
(b), 500 nm (c), and 300 nm (d) pitch reveal the increased disper- 
sivenes. The single-particle polarizability linewidth follows the local 
density of states lineshape for the vac/silver interface (e). For the 
single scatterer and two most dilute lattice cases it is feasible to com¬ 
pare the single particle linewidth to the overall lattice polarizability 
linewidth. For 1000 nm pitch, the single particle contribution is of 
comparable strength, while for 500 nm pitch, the renormalization ef¬ 
fect stemming from the lattice dominates the linewidth (f). 


response to the incident (as opposed to total) field, is renor¬ 
malized by the direct interactions with all other particles, and 
through all the interactions that take place via the mirror. 
Through the help of Fig. we examine each of these renor¬ 
malization steps. 

Figure shows the frequency and distance-dependent 
imaginary part of the dynamic polarizability for a single scat¬ 
terer and lattices with different pitch in front of the silver 
mirror. We focus on the imaginary part of the polarizabil¬ 
ity, because in free space 47r/clm(a) is associated to the ex¬ 
tinction cross section. For the bare scatterer Im(a) presents a 
Lorentzian peak at 1500 nm with peak extinction of 0.3 jive? 
corresponding to an archetypal plasmonic scatterei®. Once 
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placed in front of the mirror the single particle polarizabil¬ 
ity presents a clearly defined resonance with resonance width, 
center frequency and strength that vary with distance due to 
backaction (Fig|^ (a)). For very small separations less than 
100 nm, the resonances strongly shift when approaching the 
mirror. However, since here the dipole approximation is not 
valid we omit this range from our plots. The single-particle 
polarizability linewidth extracted from Lorentzian fits at each 
separation follows the local density of states ImG(r, r) mod¬ 
ulation associated with the vacuum/silver interface, Fig.[^(e). 

We now turn to the effective polarizability per object for lat¬ 
tices ranging from dilute to dense (1000 nm to 300 nm pitch), 
Fig.|^(b)-(d). The polarizability reduces strongly and broad¬ 
ens significantly with increasing density when considering lat¬ 
tices in absence of the mirror (result not shown). This broad¬ 
ening with i ncreasing density is well known for metamaterial 
latticeJ^^I^, and is due to the fact that many dipoles oscillating 
in phase (A:|| = 0) experience superradiant damping (radiated 
power of N dipoles being proportional to N^). Another view¬ 
point is that the extinction cross section Anklm^a) per object 
must remain below the unit cell area. Once placed in front of 
the mirror, we note that when going from dilute to dense lat¬ 
tices the polarizability is increasingly dominated by a strong 
feature that is neither solely attributable to the interface nor 
to the lattice. The weak fringes in polarizability for the sin¬ 
gle particle sharpen considerably, and become very strongly 
dispersive, forming features along lines of constant d/X. The 
oscillation period reduces to 750 nm, i.e., exactly half the 
single-particle resonance wavelength. For the densest lattices 
the fringes essentially span over the entire frequency axis. 
This notion is consistent with the fact that for dense bare lat¬ 
tices, the single particle polarizability looses relevance, as the 
polarizability is increasingly determined by just the unit cell 
area. For such dense lattices, the calculated effective polar¬ 
izability in effect implies that a strong dipole can be induced 
only whenever the lattice and the mirror satisfy a Fabry-Perot 
resonance condition. For these conditions, the areal density 
of maximum achievable polarization is high, providing a clear 
indication that conditions of perfect absorption can be associ¬ 
ated with a strong increase in particle polarizability. 

The bottom row of Fig.j^shows the FWHM distance depen¬ 
dence for the lattice and interface normalized polarizability. 
The single-particle polarizability renormalized for the pres¬ 
ence of the interface from Fig. (e) serves as a reference in 
Fig.|^(f). For dilute lattices, 1000 nm pitch, Fig.|^(f) shows 
that the interface and lattice contribution to the linewidth are 
comparable. For denser lattices the renormalization effects 
due to the interaction of the lattice of dipoles with its own 
mirror image is very large, resulting in a modulation of the 
polarizability resonance bandwidth by a factor 5 (as opposed 
to the 25% effect expected for just a single scatterer at the 
same mirror). In contrast to the case of a single particle, this 
modification has a period exactly A/2, and a modulation depth 
that stays much stronger for larger distances, indicative for the 
fact that backaction is now through an entire extended lattice 
of mirror dipoles. Importantly, while a single particle mir¬ 
ror image yields backaction through a spherical wave, for a 
mirror-image of scatterers forming a periodic array, backac¬ 


tion is through a reflected wave at k\\ =0 stemming from com¬ 
bined and fixed-phase related outgoing and receiving refiected 
plane waves. 

An apparent paradox is that the renormalization of the lat¬ 
tice polarizability through backaction with its own mirror im¬ 
age is very large, while the resulting renormalization of the 
effective lattice refiection constant is small. Qualitatively, the 
renormalization results in the highest polarizability right when 
the lattice and its mirror image form a subradiant pair as ra¬ 
diative damping is weakest. This condition is exactly coinci¬ 
dent with the incident wave Eq having a null, meaning that 
the reflectivity is least sensitive to the lattice right when the 
polarizability is most strongly renormalized. From this rea¬ 
soning, it is evident that refiection contrast is a poor probe of 
the large variations in polarizability. While back action effects 
observed in the far field measure of reflectivity are weak, we 
expect near field measures such as absorption, nonlinear fre¬ 
quency conversion, or fiuorescence spectroscopies to feel the 
polarizability enhancement much more strongly. 


V. CONCLUSION 

To conclude, we have derived a calculation method to pre¬ 
dict the optical response of metasurfaces and particle gratings 
held in front of a silver interface. We have applied this model 
to study the backaction of the interface on the polarizability of 
the lattice. We conclude that while a Fabry-Perot which takes 
as input the reflectivities of the substrate and lattice alone pro¬ 
vides a reasonably satisfactory description of the joint optical 
response, in fact there is a surprisingly large effect of the sub¬ 
strate on the particle polarizability. This could significantly 
impact the result of such composite systems in the near field. 

Finally, we note that the calculation method is easily ex¬ 
tended to lattices embedded in any arbitrary stratified system. 
Indeed, for any stratified system the Green function is analyt¬ 
ically known and always takes the form of the expression in 
square brackets in Eq. m albeit with different tensorial pref¬ 
actor M. However, the key point is that the difficult paral¬ 
lel wave vector integral always reduces to a discrete sum that 
is easy to evaluate. Thereby, we expect that this model will 
be of large utility not just for metasurface physics, but also 
for predicting the physics of diffractive outcoupling structures 
patterned on and in LEDs and remote phosphors to improve 
light generation and directional extraction. Importantly, our 
model is not restricted to non-diffractive systems, nor to any 
particular incident k \\. Thereby, it contains the full richness of 
grating diffraction, distributed feedback, and hybridization of 
lattice modes with the dispersion of any guided modes in the 
stratified system. 
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Appendix A: Illustrative experiment 

To illustrate the problem described in this manuscript, we 
present a simple experiment. We fabricated squared periodic 
arrays of Au scatterers using electron beam lithography in 
ZEP-520 resist, thermal evaporation of Au, and lift-off. The 
arrays were covered by a MicroposiE^ S1813G2 UV resist 
wedge by gray-scale lithography, which was used by us earlier 
as dielectric spacer of smoothly varying height for Drexhage 
experiments^. The lattice is located in a quasi-homogeneous 
environment at a glass (n = 1.52) - resist (n = 1.6 interface. 
The scatterers are Au split rings on glass with 100 nm arm 
length, 50 X 57 nm gap and 30 nm height as measured by elec¬ 
tron microscopy. While our experiment was originally con¬ 
ceived to probe magnetic back action effectJ^, according to 
our polarizability tensor calculation^^ such small split rings 
scatter strongly as electric dipoles (polarizable only along the 
gap, resonant around 1300 nm), with a negligible magnetic 
response. Thus they are essentially equivalent to plasmon 
nanorods. The entire structure was coated with 100 nm sil¬ 
ver that acts as mirror. The sample structure is presented 
in Fig. (a) together with a bright-held low-magnihcation 
microscopy image Fig. (b). In this image, the silver sur¬ 
face appears bright white, while periodic arrays of SRR lat¬ 
tices with 150 nm periodicity are visible as darker squares 
(50 /im e-beam write helds) arranged in a checkerboard pat¬ 
tern. The millimeter-sized wedge was printed over a macro¬ 
array of split ring arrays, as shown in Fig.[^(b). The wedge is 
faintly visible, it starts at the lower end of the image (indicated 
by marker) and slopes upwards along the vertical coordinate 
(x-direction) in the image, which results in a a smooth color 
change of the arrays. 

We performed normal incidence rehectivity measurements 
with a weakly focused beam from a halogen lamp, polar¬ 
ized along X to match the split ring fundamental resonance. 
Given the 5 /im spot size, we can measure as function of 
spacer thickness by taking several measurements in a write- 
held, and furthermore by moving from held to held along 
the wedge. White light is spectrally resolved with a hber 
coupled NIR-spectrometer (Acton SP-2350i, coupled to an 
OMA V 1024 pixel InGaAs array). We plot the measured 
rehectivity for a very dense (150 nm pitch), and a more dilute 
(300 nm pitch) lattice. Fig. |^(c,d). The bare split ring reso¬ 
nance (no mirror, split ring embedded in the dielectric) is at 
around cc = 1.5 x 10^^ (1250 nm), which is consistent with 

reports by de Hoogh et al. in Ref. |4T] For the dilute lattice we 
retrieve a high overall rehectivity, with a set of very deep min¬ 
ima around the single particle resonance. For the dense lattice, 
the rehectivity shows strong contrast across the entire spectral 
bandwidth and clearly asymmetric fringes that hip asymmetry 
at the single particle resonance. 

The standout features of interest for this paper are the fol¬ 
lowing. First, extremely large rehectivity contrast is obtained 


3 sample sketch 




C 150 nm pitch 



d 300 nm pitch 



0.2 0.4 0.6 0.8 1.0 
distance (pm) 


FIG. 5. The sample sketch (a) depicts the following sample structure: 
lattices consisting of gold split ring resonators in a periodic arrange¬ 
ment on glass are separated from a silver mirror by a dielectric wedge 
(Microposit™ S1813G2 UV-resist) (not to scale). Under bright held 
illumination (b) these lattices feature a smooth color change in x- 
direction revealing the smooth change in lattice-mirror separation on 
a /xm-length scale. (c,d) Experimental rehectivity data for split ring 
lattices with a pitch of (c) 150 nm and (d) 300 nm as a function of 
angular frequency and lattice-mirror separation. 


with both lattices. Although the rehectivity of the bare lattice 
is on the order of a few percent (not shown), with about 30% 
single pass absorption, the metasurface etalons show both near 
zero-percent and near-perfect rehection behavior for particu¬ 
lar distances. Second, the two different densities give remark¬ 
ably different responses. For the more dilute case, the res¬ 
onance yields pockets of deeply suppressed rehectivity, i.e., 
strong absorption coincident with the resonance. In contrast 
in the denser case, the rehectivity is high and featureless right 
at split ring resonance, but shows deep troughs either just to 
the red or just to the blue of a geometrical resonance con¬ 
dition (thickness an integer times a half-wavelength, mean¬ 
ing that the array is in a node of the standing wave that the 
incident wave forms at the mirror). The asymmetry hips as 
one goes through the particle resonance. A somewhat similar 
phenomenon as we observe for the dense lattice was reported 
by Ameling et al. for a lattice in a microcavity, who inter¬ 
preted their observation as an anticrossing of a Fabry -Perot 
resonance and a resonance of the scattering objects^^^. 


reflectivity reflectivity 
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Appendix B: Green’s function 


2. Lattice-summed Green function 


1. Green’s function near a planar interface 

As described by Eq. 0 in Section |I^ the dyadic Green’s 
function near a planar interface can be separated into a free 
space (in absence of an interface) and a reflected part 

G(r,r') = Gfree(Ar') + Grefl(r,r'). 

We start from the free space scalar Green function 

gi/c|r-r'| 


G(r,r') = 


r - r' 


(Bl) 


where k = ujujc is the wave number in the medium of index 
n that contains source r' and observation point r. One finds 
dyadic Green functions by differentation. The 3x3 dyadic for 
free-space is given by 


Gfree(r, r') = + Vv] G(r, r') 


(B2) 


One finds the reflected Green function in three steps. First, 
one rewrites Eq. through the angular spectrum represen¬ 
tation that states 


Jk\r\ 


— [ dnTpiq-rilpife.kl 

2% J ^ 


where q represents in-plane wave vector, kz = \/k?‘ - |qp and 
r = (r||, ^). As second step one applies the differentiation op¬ 
erator ( |B2| ), while as third step one projects the result in s- and 
p-polarization, and applies the Fresnel reflection coefficient. 
The result reads 




ikz\z+z'\ 


(qI + Qv) 


M,,- 




(B3) 


with rs,p of course dependent on q and 




% -Qxqy ot 
^ 0 


Qx Qy 
0 


Qx 

0 


0 


and 


Mp = 


Qlk'^z Qx^ykl qxkzH^^ 
qxQykl qykl qykzhl'^ 
-qxkz\q\‘^ -qykz\q\‘^ j 


The free space lattice sum is the summation of the dyadic 
Green function over all real-space lattice sites Rmn 


= E G(R„„ -r||). 


^Zkii -Rrr 


which is poorly converging. Following Ref. |42|one can split 
this sum into two convergent terms and following 
a technique pioneered by Ewald. Here we describe the imple¬ 
mentation for summing the scalar Green function. Consider 


r(fc||,r)= ^G(R™„-r> 


Akn-Flrr 


(B4) 


with G(Rmn “ I*) = G(Rmn 51 *) given in Eq. (Bl). It is pos¬ 
sible to write 


with 


r(fe||,r)=r(i)(fc||,r) + r( 2 )(fe||,r) 


TT f )’^ll 




k^ 




and 


r( 2 )(fe||,r)=^ 


^-'ik% 


AkwRrr 


Nerfc^^-|z|7?j I 


(B5a) 


2pr} 




+ e 


-ikprr 


'erfc(pmnr/-^) } 


(B5b) 


r = 


Here we used k = cj/c, 

and pmn = |Rmn “ '^ 11 1- Note that the first summa¬ 
tion is in reciprocal space, while the second is in real space. 
The convergence of Eq. ( |B5b| ) and Eq. ( |B5a| ) can be opti¬ 
mized by setting the parameter 77 around p = \/7Tla, where 
a is the lattice constant. The limit of the summation over m 
and n should be at least bigger than the expected number of 
propagating grating diffraction orders. As discussed in Refill, 
in our metamaterials calculations, with essentially no grating 
orders, i.e., ka < 27r, we already obtained converged lattice 
sums for |m, n| <5. The extension to the required dyadic lat¬ 
tice sum directly follows by applying the differential operator 
in ( |B2l l. Helpful diferentiation rules that exploit the spheri¬ 
cal resp. cylindrical coordinate formulation of ( B5b|B5a ) are 
listed in|23l 
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